BAR4R-CONF-04/26 
SLAC-PUB-10611 
July 2004 



A Search for r^Ur Recoiling Against B D*^(. vn 

The BABAR Collaboration 
February 7, 2008 



Abstract 

We present a search for the decay — > t^v^- in 124.1 X 10^ r(45) decays recorded with 
the BABAR detector at the SLAC PEP-II i?-Factory. A sample of events with one reconstructed 
exclusive semi-leptonic B decay {B^ D*^t~V^) is selected, and in the recoil a search for B~^ — > 
T^Ur signal is performed. The r is identified in the following channels: t+ e^v^^^ — s- ^'^I'^Vr, 
r+ t^'^Vt, r"*" — > -K^-K^Vr, r"*" — > 7r"'"7r~7r"^I7^. We find no evidence of signal, and we set a 
preliminary upper limit on the branching fraction of B{B~^ t^^t) < 4.3 x 10~^ at the 90% 
confidence level (CL). This result is then combined with a statistically independent BABAR search 
for B^ — > T^u-j- to give a combined preliminary limit of B{B~^ — > t^Ut) < 3.3 x 10~^ at 90% CL. 
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1 INTRODUCTION 



In the Standard Model (SM), the purely leptonic decay t~^Ut ^ proceeds via quark annihi- 

lation into a boson (Fig. Its amplitude is thus proportional to the product of the i?-decay 
constant Jb and the quark-mixing- matrix element Vub- The branching fraction is given by: 



T V) 



G2 2 

Stt 



1 



m 



B 



tb+Ib 



Vub 



(1) 



where we have set fi = c = 1, Gp \s the Fermi constant, V^b is a quark mixing matrix element ^H], 
/b is the meson decay constant which describes the overlap of the quark wave-functions inside 
the meson, is the lifetime, and and rrij- are the B^ meson and r masses. This expression 
is entirely analogous to that for pion decay. Physics beyond the SM, such as supersymmetry or 
two-Higgs doublet models, could enhance B{B^ — > t'^Vt) by up to a factor of five through the 
introduction of a charged Higgs boson [S;. 

Current theoretical values for (obtained from lattice QCD calculations) have large uncer- 
tainty, and purely leptonic decays of the B~^ meson may be the only clean experimental method of 
measuring precisely. Given measurements of \ Vub\ from semi- leptonic processes such as B ^ Triu, 
fs could be extracted from the measurement of the B^ t~^Vt- branching fraction. In addition, 
by combining the branching fraction measurement with results from B mixing, the ratio |Vub|/|Vt(i| 
can be extracted from B{B'^ — > t^v.,-') j where Am is the mass difference between the heavy 
and light neutral B meson states. 




Figure 1: Purely leptonic B decay B^ — > t^v^- proceeds via quark annihilation into a boson. 

The decay amplitude is proportional to the lepton mass and decay to the lighter leptons is 
severely suppressed. This mode is therefore the most promising for discovery at existing experi- 
ments. However, challenges such as the large missing momentum from several neutrinos make the 
signature for B^ — > t^v^- less distinctive than for other leptonic modes. 

The SM estimate of this branching fraction is (9.3±3.9) x 10"^ using \Vub\ = (3.67±0.47) x 10"^ 
and /s = 0.196 ± 0.032 jH in Eq. d 

Purely leptonic B decays have not yet been observed. CLEO 5 and experiments at LEP [HIITIIHI 
have searched for this process and set limits on the branching fraction at the 90% CL. The most 

^Charge-conjugate modes are implied throughout this paper. The signal B will always be denoted as a decay 
while the semi-leptonic B will be denoted as a B~ to avoid confusion. 
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stringent limit on B{B^ t^i^t) comes from the L3 experiment jS]: 

]3{B+ T+Ur) < 5.7 X 10"^ at the 90% CL. (2) 

The BABAR collaboration has results on the search for B^ — > r^z^r decays [Oj using a sample of 
88.9 X 10^ T{4S) decays; the upper limit is: 

B{B+ T+Vr) < 4.2 X 10"^ at the 90% CL. (3) 

In this paper we report on a preliminary result from a different analysis on a larger dataset. 

2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II storage 
ring. The sample corresponds to an integrated luminosity of 112.5 fb~^ at the T(4S) resonance 
(on-resonance) and 11.9 fb~^ taken 40MeV below the T(4S) resonance (off-resonance). The on- 
resonance sample consists of about (124.1 it 1.4) x 10^ T(4S') decays (BB pairs). The collider is 
operated with asymmetric beam energies, producing a boost of (3^ ~ 0.56 of the T(45) along the 
collision axis. 

The BABAR detector is optimized for asymmetric energy collisions at a center-of-mass (CM) 
energy corresponding to the T(45) resonance. The detector is described in detail in reference |10| . 
Charged particle tracking is provided by a five-layer double-sided silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH) contained within the magnetic field of a 1.5 T superconducting 
solenoid. The tracking system provides momentum reconstruction of charged particles and measures 
energy loss {dE/dx) for particle identification. Additional charged K~tt particle identification is 
provided by a ring- imaging Cherenkov detector (DIRC), which exploits the total internal reflection 
of Cherenkov photons within the purified quartz bars. The energies of neutral particles are measured 
by an electromagnetic calorimeter (EMC) composed of 6580 CsI(Tl) crystals, which provides an 
energy resolution aE/E = {2.2,/ E^/"^ 1.9)%, where E is in GeV. For 20 MeV and 1 GeV clusters 
the EMC energy resolutions are ~6% and ~3%, respectively. The magnetic flux return of the 
solenoid (IFR) is instrumented with resistive plate chambers in order to provide muon and neutral 
hadron identification. 

A GEANT4-based Monte Carlo (MC) simulation is used to model the signal efficiency 
and the physics backgrounds. Simulation samples equivalent to approximately three times the 
accumulated data were used to model BB events, and samples equivalent to approximately 1.5 
times the accumulated data were used to model e+e~ — > uu, dd, ss, cc, and r+r~ events. A large 
sample of signal events is simulated, where one of the B meson decays to B"*" — > t^Vt- Beam 
related background and detector noise from data are overlaid on the simulated events. 

3 ANALYSIS METHOD 

Due to the presence of multiple neutrinos, the B+ t^v^- decay mode lacks the kinematic con- 
straints which are usually exploited in B decay searches in order to reject both continuum and BB 
backgrounds. The strategy adopted for this analysis is to reconstruct exclusively the decay of one 
of the B mesons in the event, referred to as "tag" B, and to compare the remaining particle(s) in 
the event, referred as the "signal side", with the signature expected for the decay B"*" — > t~^Ut. In 
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order to avoid experimenter bias, the signal region in data is not examined ("blinded") until the 
selection is optimized based on MC simulation. 

The tag B is reconstructed in the set of decay modes D*^i~Vi, where ^ is e or fi. The D*^ 

is reconstructed in D^tt^ and D^j modes. The is reconstructed in four decay modes: K~7r~^, 
K~TT~^TT~TT~^ , K~TT~^TT^, and K^TT~^TT~. The is reconstructed only in the mode — > tt~^tt~ . On 
the signal side the B~^ — > t'^Vj- signal is searched for in both leptonic and hadronic r decay modes: 
e'^Vei^T-, fJ'~^i^fjpTi T^^T^Ti ^ tt^tt^Vt, TT^ TT^ tt^Vt-. The branching 

fractions of the above r decay modes are listed in Tabled Most of the kinematic variables used 
for event selection and background rejection are measured in the CM frame. 

Table 1: Branching fractions for the r decay modes used in the i?"*" — > t'^i't- search 



Decay Mode 


Branching Fraction (%) 




17.84 ± 0.06 


T+ M+l/^Fr 


17.36 ± 0.06 




11.06 ± 0.11 




25.42 ± 0.14 




9.16 ± 0.10 



3.1 TAG B RECONSTRUCTION 

The tag B reconstruction proceeds as follows. First we reconstruct the D*^ candidates in the above 
four decay modes using tracks and/or a vr*^. The tracks are required to meet particle identification 
criteria consistent with the particle hypothesis, and are required to converge at a common vertex. 
The TT^ candidate is required to have invariant mass between 0.115-0.150 GeV/c^ and its daughter 
photon candidates must have a minimum energy of 30 MeV. The mass of the reconstructed 
candidates in K~7r~^, i('~7r^7r~7r"^, and Kg7r~^7r~ modes are required to be within 40 MeV of the 
nominal mass In the K^tt^tt^ decay mode the mass is required to be within 70 MeV of the 
nominal mass (4,. 

The D*^ candidates are reconstructed by combining the candidates with a soft vr*' or 7, 
whose momentum in the CM frame is less than 0.45 GeV/c. The mass difference between D*^ and 
D° (AM) is restricted to be within 0.13-0.17 GeV/c^ and 0.12-0.17 GeV/c^ for D°7r° and D^j 
modes, respectively. We further require that the photon used in D*^ — > D^'y reconstruction has a 
minimum energy of 100 MeV. 

Finally D*^£ candidates are reconstructed by combining D*^ with an identified electron or muon 
of a momentum above 1.0 GeV/c in the CM frame. The D*^ and i candidates are required to meet 
at a common vertex. An additional kinematic constraint is imposed on the reconstructed D*^£ 
candidates: Assuming that the massless neutrino is the only missing particle, we calculate the 
cosine of the angle between the D*^£ candidate and the B meson, 

. 2EBEj^,oi -ml- m^,o^ 
cos6Ib_^.o^ = . (4) 

Here {E£)*o£, Pd*o£) and {Eb, Pb) are the four-momenta in the CM frame, and m£)»o^ and niB are 
the masses of the D*^£ candidate and B meson, respectively. Eb and the magnitude of pB are 
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calculated from the beam energy: Eb = -Ebeam = ^cm/2 and \pb\ = y-E'beam ~ "^s- Correctly 
reconstructed candidates populate the range [-1,1], whereas combinatorial backgrounds can take 
unphysical values outside this range. We retain events in the interval —1.1 < cos6'^_£,*0£ < 1.1, to 
take into account the detector energy and momentum resolution. 

If more than one suitable D*^i candidate is reconstructed in an event, the best candidate is 
selected based on the mass and AM. A two dimensional likelihood function is formed by taking 
the product of the mass and AM distributions obtained from MC simulation. We select the 
candidate with the largest likelihood. 

The following additional cuts are applied on the selected best candidate. The AM value for 
the best candidate is required to be between 0.135-0.150 GeV/c^ and 0.130-0.155 GeV/c^ for 
candidates reconstructed in D^tt^ and D^^ modes, respectively. The angle between the and 
the soft vr'' or 7 from D*^ decay in the CM frame is restricted to be less than 60° and 90° for D*^ 
reconstructed in D^n^ and D^^ modes, respectively. The sum of the charge of all the particles in 
the event (net charge) must be equal to zero for events with selected candidates. 

At this stage of the selection, the observed yield in data and the predicted yield in the MC 
simulation agree to within approximately 6%. This discrepancy is corrected by scaling the yield 
and efficiency obtained from MC simulation. The scale factor of 0.937 is used to correct the tag 
B reconstruction efficiency in the signal MC simulation. The systematic error associated with this 
correction will be described in Sec. |SJ The corrected tag reconstruction efficiency in the signal 
MC simulation is (1.818 ± 0.074) xlO^^. Fi gure[2l shows the AM distributions of the selected best 
D*^ev candidates. 



3.2 SELECTION OF B+ t+v^ DECAYS 

After the tag B reconstruction, in the signal side the r from the B^ t^Ut decay is identified 
in one of the following modes: r"*" — > e^Vei'Ti fJ'^i^fjpT, t^^T^t, tt~^tt^Vt-, or 

TT^TT^ TT^u-T-. We select events with one or three signal-side track(s). The event is rejected 
if any of the signal-side tracks fail the following selection criteria: it must have at least 12 DCH 
hits, its momentum transverse to the beam axis, pT, is greater than 0.1 GeV/c, and its point of 
closest approach to the interaction point is less than 10.0 cm along the beam axis and less than 



1.5 cm transverse to the beam axis. Figure 3(a) shows the distribution of the number of signal- 
side tracks before this cut. The invariant mass of a signal-side vr*^ candidate must be between 
0.10-0.16 GeV/c^, the shower shape of the daughter photon candidates must be consistent with 
an electromagnetic shower shape and the photons must have a minimum energy of 50 MeV. 

The most powerful variable for separating signal and background is the remaining neutral energy 
(£^extra). Calculated by adding the CM energy of the photons that are not associated with either the 
tag B or the vr*^ candidate from r"*" tt'^tt^Vt signal decay. The photon candidates contributing to 
the -Eextra Variable have minimum cluster energies of 20 MeV. For signal events the neutral clusters 
contributing to -Ecxtra can only come from processes like beam-background, hadronic split-offs, and 
bremsstrahlung. Therefore the signal events peak at low -Eoxtra values and the background events, 
which contain additional sources of neutral clusters, are distributed towards higher -Eextra values 
(see fig. |3(b)] ). The E^extra < 0.3 GeV region is defined as the signal region. This -Eextra < 0.35 
GeV region, which is slightly larger than the signal region, is kept blinded in on-resonance data 
until the selection is optimized. 

The different signal modes are distinguished by the following selection criteria. The r+ 
e^VgUr, r"*" — > ix^u^Vt, — 'k'^Vt, and r"*" tt'^'k^Vt- signal modes, all of which contain one 
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(a) The AM distribution of D*°{D°-K°)ev can- 
didate plotted for — » t^Vt signal simulation 
(top) and for data and background simulation 
(bottom). 



(b) The AM distribution of D*°{D°'y)ev candi- 
date plotted for —* t^Vt signal simulation 
(top) and for data and background simulation 
(bottom). 



Figure 2: The AM distributions of D*^ev candidate with D*^ reconstructed in D^tt^ (fig. 2(a) I 
and D^j (fig. 2(b) ) modes, plotted for T~^Ur simulation, background simulation and data. 

The AM, and the angle between the and the soft vr*^ or 7 requirements are not applied on 
the 0*^61^ candidates plotted in these distributions. Here the simulated background events are 
scaled to match the data yield. The normalization of the signal MC events is arbitrary. These 
distributions are similar for the D*^^v candidates. 



charged track, are separated by particle identification. Both the r"*" — > tt+z/t- and the — > vr+vr 17,- 
modes contain a pion signal track and are characterized by the number of signal-side tt^ mesons. 
The r+ vr+Tr^vr+I/.^ mode contains three signal-side tracks. These signal selection requirements 
are as follows. 

• Particle identification: 

— For the r"*" — > e^^UePr selection the track must be identified as an electron and not 
identified as a muon or a kaon. 

— For the ii^u^yj- selection the track must be identified as a muon and not identified 
as an electron or a kaon. 

— For the ir^Vr and r"*" ■K'^ir^T'r selection we require that the track is not identified 
as an electron or a muon or a kaon. 

— For the tt^tt^tt^Vt selection each of the tracks must be identified as a pion and 
not identified as an electron or a muon or a kaon. 

• Signal-side vr" multiplicity: 
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(a) Number of signal-side tracks for D*° {D''^'y)eu 
tag B candidate plotted for B'^ —^ t^Vt sig- 
nal MC simulation (top) and for data and back- 
ground MC simulation (bottom). 



(b) -Bcxtra distribution after the D*^ {D^-y)ev tag 
B candidate selection plotted for B^ — > t'^Vt 
signal MC simulation (top) and for data and 
background MC simulation (bottom). 



Figure 3: The number of signal-side tracks and -Eextra distributions after tag B selection, plotted 
for — > t~^Vt- MC simulation, background MC simulation, and data. The net charge requirement 
is not applied on the events plotted here. Here the simulated background events are normalized 
to on-resonance data luminosity. The normalization of the signal MC events is arbitrary. These 
distributions are similar for other tag B candidates. 

— For the t+ — > tt+Ft- selection we require the event to contain no signal-side vr*^. 

— For the r+ ■k'^tt^Vt selection we require that the event contains at least one signal-side 



• -Eextra requirement: 

— For all the signal modes -Eextra must be less than 0.3 GeV. 

Background consists primarily of B^B^ events in which the tag B meson has been correctly 
reconstructed. The recoil side contains one or three track(s) and the additional particles which are 
not reconstructed by the tracking detectors or calorimeters. Typically these events contain one or 
more and/or neutrinos, and frequently also additional charged or neutral particles which pass 
outside of the tracking and calorimeter acceptance. Background events also contain B^B^ events. 
The continuum background contributes to hadronic r decay modes. In addition some excess events 
in data, most likely from two-photon processes which are not modeled in MC simulation, are also 
seen. These backgrounds can be suppressed by the following constraints on the kinematics of the 
B^ — > T^u-r- decays. 



13 



• Missing mass: The missing mass is calculated as follows. 



Mmiss = ^ {Er(iS) - EvmY - (Pr(45) - PVisf- (5) 

Here {Ey[4s)^ Pr{4S)) is the four-momenta of the T{4S), known from the beam energies. The 
quantities £^vis and pVis are the total visible energy and momentum of the event which are 
calculated by adding the energy and momentum, respectively, of all the reconstructed charged 
tracks and photons in the event. 

— For the e^z/eFr and t+ ^^v^Vj- selections, events with missing mass less than 
4 GeV/c^ are rejected. 

— For the vr+F,- and r+ tt^tt^Vt selections, the missing mass is required to be 
greater than 3 GeV/c^. 

— For the r"^ ■k^tt^tt^Vt selection, the missing mass is required to be greater than 2 
GeV/c2. 

• Maximum CM momentum of the r daughter: 

The following maximum CM momentum requirements are applied to the r daughter particles. 

— The electron candidate from the r"*" — > e^VeT^r decay must have a CM momentum of less 
than 1.4 GeV/c. The CM momentum requirement is not applied to the r"*" — > fi'^VpUr 
selection because of the following reason: The momentum spectrum of the lepton from 
T decays peaks below 1 GeV/c. The particle identification efficiency for low momentum 
muons is lower than that for low momentum electrons. Therefore, applying the maximum 
momentum cut reduces the selection efficiency of the r"*" — > ^'^u^Vr mode significantly. 

— For the three hadronic r decay modes, the vr from r"^ t^~^i^t, the tttt^ combination 
from r"*" — > Tr'^TT^Ur, or the Svr combination from — > tt^ tt^ tt^Vt must all have CM 
momenta less than 2.7 GeV/c. 

The —>■ -k^tt^Vt- and r"^ -k^tt^tt^Vt decays proceed via intermediate resonances. For 
these modes further background rejection can be achieved by applying the following requirements 
on the intermediate mesons. 

• selection: 

The t'^ — s- TT^T^^v-r- decay proceeds via an intermediate p"^ state. The signal-side track is 
combined with a signal-side vr*^ to form the p"*" candidate. In events with more than one 
signal-side tt*^, the candidate with invariant mass closest to the nominal vr*^ mass [1] is chosen. 
The invariant mass of the reconstructed p'^ is required to be within 0.55-1.00 GeV/c^. A 
quantity similar to cos^^.^j.o^, which is defined in section ITTl can be reconstructed for 
T pv as follows: 

— — ' 

where {Er, Pr) and {Ep, pp) are the four-momenta in the CM frame, mr and mp are the 
masses of the r and p candidate, respectively. The quantities \pr\ and E^- are calculated 
assuming that the r is from the t~^Vt- decay, and that the is almost at rest in the 

CM frame. Candidates outside of —1.1 < cos Or- p < 1.1 are excluded. 
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• p'^ and a\ selection: 

The T decays to three charged tracks via two intermediate resonances: afv^-, —^ 

p^iT~^ 1 and — > 7r"^7r~. The 7r~7r"'" combination with an invariant mass closest to the nominal 
p^ mass 3 is selected as the best p^ candidate. The invariant mass of the selected p^ must be 
within 0.55-1.00 GeV/c^. The CM momentum of the selected p^ candidate is required to be 
greater than 0.5 GeV/c. The invariant mass of the three signal tracks must be within 1.0-1.6 
GeV/c^. The total CM momentum of the three tracks has to be greater than 1.0 GeV/c. 
The three tracks are also required to converge to a common vertex and the candidates are 
rejected if the vertex fit probability is less than 0.1%. For r aiu decay the quantity, 

= — wKi — • 

is obtained using the similar procedure used for calculating cos6't— p. Here (Sri pV) and {Eax-, 
Pai) are the four-momenta in the CM frame, m.^ and are the masses of the r and a\ 
candidate, respectively. Candidates not satisfying — 1.1 < cos^r- ai < 1-1 are excluded. 

The signal selection criteria for all five signal modes are summarized in Table 



Table 2: The selection criteria for different signal modes are listed in this table. The symbols P* 
and Mj;, used in the table, correspond to the CM momentum and invariant mass of x, respectively. 













One signal-side track 


Three signal-side tracks 


track quality requirements for each signal track 


electron 


muon 


not electron 


pion 


not muon 


not electron 


not muon 


not electron 


not kaon 


not kaon 


not kaon 


not muon 










not kaon 


none 


none 


No signal-side tt*^ 


Non-zero signal-side tt" 


none 


£;extra < 0.3 GeV 


Missing Mass > 4 GeV 


Missing 


Mass > 3 GeV 


Missing Mass > 2 GeV 


P* of 


none 


P* of 


P* of 


P* of 


signal-side e 




signal-side n 


signal-side 


of signal-side 


track < 1.4 GeV 




track < 2.7 GeV 


■Kn° < 2.7 GeV 


3 TT < 2.7 GeV 


none 


none 


none 


p"^ selection: 
0.55 < Mp+ < 1 GeV 

-1.1 < COs9r-p < 1.1 


selection: 
0.55 < A'4+^- < 1.0 GeV 
P;+^- > 0.5 GeV 
1.0 < Mz^ < 1.6 GeV 
PJ^ > 1.0 GeV 
Vertex prob. of 3 tracks > 1% 
-1.1 < COsSr-ai < 1.1 



3.2.1 SIGNAL EFFICIENCY 

The "signal-side selection efficiencies" for the r decay modes are determined from signal MC sim- 
ulation and summarized in Table 01 For each r selection mode, the signal-side efficiency (ej, where 
i = selection mode) is computed as the ratio of the number of events surviving the requirements 
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of that selection mode to the number of events where a tag B meson is reconstructed. In the 
computation of the total signal-side efficiency for each selection we take into account the cross-feed 
from other r decay modes reported in Table 01 

The selection efficiency for r"*" fi^v^r is low compared to that of the r+ — > e^fe^^r mode, 
because of the fact that the momentum spectrum of the signal muons peaks below 1 GeV/c, where 
the muon detection efficiency is low. Since no minimum momentum requirement and no tight pion 
identification criteria are applied to the r"*" — > tt^Vt signal selection, electron and muon signal tracks 
that fail particle identification requirement get selected in this mode. Any true tt'^tt^Vt signal 

events, with a missed vr'^ also get included in r"*" — > tt'^Vt selection mode. Therefore the r"*" — > tt'^T't- 
selection mode has the highest signal efficiency. Prom MC estimation we expect ~10 signal events 
at 112.5 fb-\ assuming B{B+ r+i/^) = 10~^. 

Table 3: Efficiency of the different selections (columns) for the most abundant r decay modes 
(rows). The last two rows show the total efficiency for each selection weighted by the decay 
branching fractions, and the total efficiency. The errors are statistical only. The total efficiency 
for each selection is = E^^^elfj, where el is the efficiency of the selection i for the simulated r 
decay mode j. The index j corresponds to the different r decay mode in the MC simulation, and 
fj = B{t j) are the r branching fractions from Ref. (4j. If the efficiency is zero a 90% upper 
limit is quoted. 



r decay 


eui' 










mode in 


selection 


selection 


selection 


selection 


selection 


Monte Carlo 


efficiency 


efficiency 


efficiency 


efficiency 


efficiency 


simulation 


(%) 


(%) 


(%) 


(%) 


(%) 


euV 


48.6 ± 2.8 


0.1 ± 0.1 


11.9 ± 1.2 


0.1 ± 0.1 


(< 0.2) 




0.1 ± 0.1 


25.8 ± 1.9 


53.3 ± 3.0 


0.7 ± 0.3 


(< 0.3) 


TTU 


(< 0.4) 


0.5 ± 0.3 


57.5 ± 4.0 


2.6 ± 0.7 


(< 0.4) 




(< 0.2) 


0.3 ± 0.1 


12.1 ± 1.0 


8.8 ± 0.8 


(< 0.2) 




(< 0.5) 


(< 0.5) 


1.1 ± 0.6 


(< 0.5) 


27.6 ± 3.2 




(< 0.5) 


(< 0.5) 


2.8 ± 0.8 


2.8 ± 0.8 


(< 0.5) 




(< 1.0) 


(< 1.0) 


0.5 ± 0.5 


(< 1.0) 


4.1 ± 1.5 


other 


1.2 ± 0.5 


2.1 ± 0.7 


10.0 ± 1.5 


2.5 ± 0.7 


1.0 ± 0.4 


all r decay 


8.4 ± 0.4 


4.8 ± 0.3 


21.9 ± 0.7 


3.2 ± 0.2 


2.1 ± 0.2 


Total 


40.4 ± 0.9 



3.2.2 EXPECTED BACKGROUND FROM MONTE CARLO SIMULATION 

To obtain the background estimation from the MC simulation, BB and e^e^ uu, dd, s's, 
cc, and r+r" events are scaled to equivalent luminosity in data. The estimated background in 
different selection modes is listed in Table 0] The three modes r"^ — > e'^VeUr, t"*" fi'^v^Ur, 
r"*" — > tt~^Vt- are relatively clean compared to the other two selection modes. Signal to background 
ratios in r"*" — s- e~^Uel7T-, r"*" /x"*" 1/^17,-, and r"*" — > tt'^Vj- selection modes are 10:59, 10:75, and 10:85 
respectively, assuming B{B~^ — > r+z^,-) = 10~^. For the r"*" — > tt^tt^Ut and — > tt^ tt^ ir^Uj- modes 
the signal to background ratio is 10:354 and 10:513, respectively. 
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Table 4: Expected final raw background and signal yield at 112.5 fb^ estimated from background 
and signal MC simulation. No systematic correction is applied on simulated events. The listed 
errors are statistical only. If the efficiency is zero a 90% upper limit on the expected background is 
quoted. 















B+B- 


8.29 ± 1.69 


7.25 ± 1.58 


24.18 ± 2.89 


14.16 ± 2.21 


14.16 ± 2.21 


Bm^' 


3.36 ± 0.97 


1.40 ± 0.63 


8.11 ± 1.51 


3.64 ± 1.01 


7.83 ± 1.48 


Combined 
BB 


11.65 ± 1.95 


8.65 ± 1.70 


32.29 ± 3.26 


17.80 ± 2.43 


21.99 ± 2.66 


cc 


(< 2.2) 


(< 2.2) 


9.54 ± 3.02 


8.59 ± 2.86 


2.86 ± 1.65 


uu, dd, ss 


(< 1.9) 


(< 1.9) 


(< 1.9) 


(< 1.9) 


0.81 ± 0.81 




0.48 ± 0.48 


(< 1.1) 


2.87 ± 1.17 


0.48 ± 0.48 


0(<1.1) 


Combined 
non BB 


0.48 ± 0.48 





12.41 ± 3.24 


9.07 ± 2.90 


3.67 ± 1.84 


Total Background 


12.12 ± 2.01 


8.65 ± 1.70 


44.69 ± 4.59 


26.86 ± 3.79 


25.67 ± 3.24 


Expected signal 
events for 

B{B+ T+Ur) = 

10-4 


2.03 ± 0.10 


1.16 ± 0.07 


5.26 ± 0.16 


0.76 ± 0.06 


0.50 ± 0.05 



3.2.3 BACKGROUND ESTIMATION FROM Sextra SIDE BAND IN DATA 

The -Eextra < 0.3 GeV region is defined as the "signal region" and the 0.35 < -Eoxtra < 1-0 GeV 
region is defined as the "side band" . The -Eextra shape in the MC distribution is used to extrapolate 
the data side band to the signal region. 

The number of MC events in signal region (N^^) and side band (A'^ly^g) are counted and their 
ratio (-Rmc) is obtained. 



atMC 

^MC - ^MC 
^^SideB 

Using the number of data events in the side band (N^^^^^) and the ratio -RmC; the number of 
expected background events in the signal region in data (N^^^^^) is estimated. 



Ardata Ardata d 

^^SideB - ^^SideB ' ^MC 

The background estimation for the different selection modes from the -Eextra side band extrapolation 
are shown in Table El The number of estimated background events in the signal region from the 
data side band extrapolation are in agreement with the background estimation from MC simulation 
within statistical uncertainty. 
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Table 5: Background estimation in the signal region (-Eextra < 0.3 GeV) for the different selection 
modes. 



Selection 


^^SideB 




Rmc 


Ardata 
^^SidcB 


^^sTg^ 




83.63 ± 5.67 


12.12 ± 2.01 


0.14 ± 0.03 


103.00 ± 10.15 


14.93 ± 3.05 




64.18 ± 4.83 


8.65 ± 1.70 


0.13 ± 0.03 


53.00 ± 7.28 


7.14 ± 1.80 




227.09 ± 10.86 


44.69 ± 4.59 


0.20 ± 0.02 


250.00 ± 15.81 


49.20 ± 6.39 




179.44 ± 9.37 


26.86 ± 3.79 


0.15 ± 0.02 


182.00 ± 13.49 


27.24 ± 4.57 




218.59 ± 9.55 


25.67 ± 3.24 


0.12 ± 0.02 


199.00 ± 14.11 


23.37 ± 3.53 



4 VALIDATION OF Eextra SIMULATION 

The -Eextra distribution in signal and background MC simulation are validated using various control 
samples. We compare the -Eextra distributions between on-resonance data and MC simulation, in 
both signal region and side band, using the control samples. Agreement between the distributions 
would provide validation of the -Eextra modeling in the simulation. 

4.1 £;extra IN THE SIGNAL MONTE CARLO SIMULATION 

The "double-tagged" events, for which both of the B mesons are reconstructed in tagging modes, 
B~ — > D*^i~V£ vs — > D*^£^i'i, are used as a control sample to validate the -Eextra simulation. 
Due to the large branching fraction and high tagging efficiency for these events, a sizable sample 
of such events is reconstructed in the on-resonance dataset. These double-tag events contain very 
little background due to the full reconstruction of the event. 

To select double-tag events we require that the two tag B candidates do not share any tracks 
or neutrals. If there are more than two such non-overlapping tag B candidates in the event then 
the best two are selected using the same best candidate selection criteria, as described in Sec. 13.11 
After selecting the two tag B candidates only the events with no extra charged tracks are selected. 

The -Eextra for the double-tagged sample is calculated by summing the CM energy of the photons 
which are not associated with either of the tag B candidates. The sources of neutrals contributing 
to the -Ecxtra distribution in double-tagged events are similar to those contributing to the -Eextra 
distribution in the signal MC simulation. Therefore the agreement of the -Eextra distribution between 
data and MC simulation for the double-tagged sample, in figure HJ is used as a validation of the 
-E'extra simulation in the signal MC. 

4.2 E^extra IN THE BACKGROUND MONTE CARLO SIMULATION 

The following two background control samples are used to study the agreement between data and 
simulation in the -Eextra < 0.35 GeV region. 

• Events with two remaining signal-side tracks 

• Events with non-zero net charge 
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Figure 4: The distribution of the remaining neutral energy (E'extra) for double-tagged events, plotted 
for simulation and data. 

The procedure from Sec. 13. 2. 31 the background estimation from the -Eextra side band, is used for 
this test. In TableElwe show the comparison of the number of expected data events (N^-^^^^^) in the 
signal region with the observed number of data events (N^i'g) in the signal region. The agreement 
between the above two quantities provides validation of background estimation in the low -Bextra 
region. 

Table 6: Test of background estimation in the low -Bextra region from various control samples. Signal 
selection cuts for different r decay modes are applied on the control samples. The expected and 
observed number of events in the signal region agree within error for all the samples (Pix^) = 0.15 
between the entries in the last two columns. Here = 10.77, and the number of degrees of freedom 
is 7.). 



Selection 


atMC 
^^SideB 


Ardata 
^^SideB 


Sig 




ATObs 


Two signal-side track control sample 


evV 


223.34 ± 8.68 


203.00 ± 14.25 


76.91 ± 5.12 


69.90 ± 7.29 


64.00 ± 8.00 




144.80 ± 7.14 


123.00 ± 11.09 


57.68 ± 4.43 


49.00 ± 6.29 


55.00 ± 7.42 




1132.39 ± 21.64 


1050.00 ± 32.40 


401.89 ± 12.53 


372.65 ± 17.83 


422.00 ± 20.54 




1175.99 ± 21.78 


959.00 ± 30.97 


241.36 ± 9.70 


196.83 ± 10.78 


232.00 ± 15.23 


Qnet 7^ control sample 


evV 


16.86 ± 2.33 


8.00 ± 2.83 


4.10 ± 1.14 


1.94 ± 0.92 


5.00 ± 2.24 


[ivV 


6.66 ± 1.46 


2.00 ± 1.41 


2.29 ± 0.87 


0.69 ± 0.57 







68.62 ± 5.31 


63.00 ± 7.94 


23.26 ± 2.83 


21.36 ± 4.09 


24.00 ± 4.90 




63.14 ± 5.14 


44.00 ± 6.63 


9.14 ± 1.74 


6.37 ± 1.63 


6.00 zb 2.45 
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5 STUDIES OF SYSTEMATICS 



The main sources of uncertainty in the determination of the t~^Vt- branching fraction are 

the fohowing: 

• We divide the uncertainties in signal efficiency determination into two categories. 

— Uncertainty in tagging efficiency determination 

— Uncertainty in determination of the efficiency Ei for each selection mode. 

• Uncertainty in the determination of the number of expected background events in the signal 
region for each selection mode. 

A small uncertainty of 1.1% also enters the branching ratio limit calculation from the estimation 
of the number of B^B~ events present in the data sample 

5.1 TAGGING EFFICIENCY SYSTEMATICS 

The tagging efficiency and yield in signal simulation is corrected using the "double-tagged" events. 
The selection of "double-tagged" events is described in sec. 14.11 
The number of double-tagged events (A''2) is given by 

N2 = e^N (8) 

where e is the efficiency and is the original number of events. The double tag yield in data 
is 407.0 lb 20.2. In the simulation, we find 434.4 it 12.4 double-tags, with very few non-B^ B~ 
events. The expected number of B~^B~ events in the dataset is A'^ = (62.06 it 0.68) x 10® events. 
Calculating the efficiency using the numbers of double-tags in data and normalized MC in equation 
IHl we find the efficiencies edata and eMC- The correction factor, ratio of the efficiencies between 
data and simulation, from this method is 

= 0.969 ± 0.029. 

The correction factor obtained from the double tags agrees within error with the normalization 
factor used in section IXTI to correct tag B meson yield in simulation. Therefore, the double tags 
provide the validation for the tag B yield correction. We take the 3.1% error obtained from the 
double tag method as the systematic error associated with the tagging efficiency and yield correction 
in Monte Carlo. 

5.2 UNCERTAINTIES IN THE SIGNAL SELECTION EFFICIENCIES IN 
EACH SELECTION MODE 

Besides tagging efficiency uncertainty, the contribution to the systematic uncertainties in the deter- 
mination of the efficiencies comes from systematic uncertainty on the tracking efficiency, particle 
identification, and simulation of the neutral clusters in the calorimeter which contribute to the 
-Sextra distribution. Uncertainty in the vr'' reconstruction efficiency introduces an additional contri- 
bution to the systematics in the r+ —> ir^Vr (5.0%) and r+ -k^t^^Vt (7.7 %) selection modes. 
The different contributions to the systematic uncertainty on the selection efficiencies are listed in 
table 
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Table 7: Contribution to the systematic uncertainty on the signal selection efficiencies in different 
selection modes. 



beiection 
modes 


J. 1 ■ 

trackmg 

(%) 


Farticie 
Identification 

(%) 


AT J. 1 

iNeutrai 
Reconstruction 

(%) 


iotal 
Systematic 
Error (%) 


Correction 
Factor 




1.4 


0.1 


3.1 


3.4 


0.99 




1.4 


3.0 


2.3 


4.0 


0.89 




1.4 


0.5 


3.1 e 5.0 


6.1 


1.02 




1.4 


0.1 


2.9 e 7.7 


8.3 


0.95 




4.2 


0.3 


4.7 


6.3 


1.00 



5.3 UNCERTAINTIES IN THE BACKGROUND ESTIMATION 

The background estimation is performed by extrapolating the number of events in the -Eextra side 
band in data into the signal region as described in sec. 13.2.31 The major uncertainty related 
to background estimation comes from the data and MC statistics. The modeling of the -Eextra 
variable in the background MC contributes to additional systematic corrections to the background 
estimation. The systematic corrections due to -Eextra modeling for different modes are the following: 
(1.02 ± 0.04) for r+ e+zv^IJ^, (1.13 ± 0.06) for r+ ^ n~^UfJ7r, (1-12 ± 0.03) for r+ n+Vr, 
(1.09 ± 0.04) for T+ vr+vrV and (1.07 ± 0.03) for t+ vr+vr-vr+I/^. 

6 RESULTS 

After finalizing the signal selection criteria, the signal region (-Bextra < 0.3 GeV) in the on-resonance 
data is examined. Table |S1 lists the number of observed events in on-resonance data in the signal 
region, together with the expected number of background events in the signal region. Figure 121 
shows the -Bextra distribution in data and simulation. In almost all the modes the observed number 
of events is in agreement with the expected number of background events. 

Table 8: Shown are the observed number of on-resonance data events in the signal region, together 
with number of expected background events. The background estimations include systematic cor- 
rections referred to in Sec. 15.31 



Selection 


Expected 


Observed Events 




Background Events 


in On-resonance Data 




15.15 ± 3.14 


13 




8.05 ± 2.07 


10 




55.30 ± 7.37 


72 




29.80 ± 5.10 


30 




25.10 ± 3.87 


26 



We determine the t^v^- branching fraction from the number of signal candidates Si 

in data for each r decay mode, according to si = N^-^B^B'^ — T+z^T-)£tag£i- Here N^-^ is the 
total number of BB pairs in data, Stag is the tag reconstruction efficiency in signal MC; is 
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(a) The distribution of the 
Bcxtra for events satisfying 
r+ — > e^VeT'T selection. 



(b) The distribution of the 
-Boxtra for events satisfying 
Ij.'^v^Vr selection. 



(c) The distribution of the Ecxtm 
for events satisfying tv^Vt 
selection. 




(d) The distribution of the 
Bcxtra for events satisfying 
r+ tt+tt^I't- selection. 



(e) The distribution of the -Bcxti 
for events satisfying 
7r'''7r~7r"''77^ selection. 



Figure 5: The distribution of the -Eextra for events passing the selections for different r decay 
modes, plotted for background and signal MC simulation and on-resonance data. Events plotted 
here are required to pass selections for corresponding signal modes. The background MC events 
are normalized to on-resonance data luminosity. The signal MC events for the r+ e^VePr, 
r"*" — > fi'^UfjU-r, and r"*" tt'^Vt modes are normalized assuming B{B^ — > t^Vj-) = 10"'^, whereas 
for the normalization of the signal MC events in the — > -k'^'k^Vt and r+ vr+vr^vr+FT- modes 
we use B{B^ — > t~^Vt) = 10-2. The above B{B+ — > T^Ur) assumptions for the normalization of 
the signal MC events are used to compare the shapes of the i^extra distributions between signal MC 
events, background MC simulation, and on-resonance data events. 
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the signal-side selection efficiency in different r decay modes calculated with respect to the total 
number of reconstructed tag B mesons. TableElshows the values of N^-^, etag and £i, after applying 
appropriate systematic corrections (sec. ISJ. 

Table 9: Shown are the values for the quantities -/V^-g, etag and ej, after applying systematic 
corrections. 





Value 


Total Uncertainty (%) 




(124.1 ± 1.4) X 10** 


1.1 


£tag 


(1.82 ± 0.074 ± 0.055) x 10"^ 


5.06 




(8.36 ± 0.42 ± 0.28) % 


6.03 




(4.30 ± 0.28 ± 0.17) % 


7.56 




(22.34 ± 0.72 ± 1.36) % 


6.91 


^TT+TT^T't 


(3.01 ± 0.24 ± 0.25) % 


11.45 




(2.07 ± 0.20 ± 0.13) % 


11.53 



The results from each decay mode are combined using the ratio Q = C{s + b)/C{b), where 
C{s + b) and C{b) are the likelihood functions for signal plus background and background-only 
hypotheses, respectively [TH [T3] : 

c{s + 6) ^ n - — n"^^' ^(^) - n (9) 

Since we have no evidence of signal we set an upper limit. The statistical and systematic 
uncertainties on the expected background (6j) are included in the likelihood definition by convolving 
it with a Gaussian distribution {Q). The mean of ^ is 5j, and the standard deviation ((T5. ) of Q is 
the statistical and systematic errors on bi added in quadrature [TB] . 

C{si + bi) ^ C{si + bi) g{bi, dfej (10) 

(similarly for C{bi)). We determine the limit on the branching fraction to be B{B'^ — s- r+i/^) < 
4.3 X 10^^ at the 90% C.L. Figure El shows the likelihood ratios as a function of B{B~^ T^Vr). 
The solid curve corresponds to the case in which the uncertainty on the expected background is 
not included. The effect of including the uncertainty on the expected background can be seen from 
the dashed curve. 

The measured branching fraction, the value that maximizes the likelihood ratio estimator is 
1.9l}:? X 10"^- 

The BABMi Collaboration performed also a search for the T~^Vr decay 0, where the tag B 

mesons are reconstructed in hadronic modes B~ D^*^^X~ . Here X~ represents a combination of 
up to five charged pions or kaons and up to two vr'' candidates. The hadronic reconstruction analysis 
is statistically independent from the current analysis and has obtained a limit B{B^ — >■ T^fr) < 
4.2 X 10~^ at the 90% C.L. We combine the results from the statistically independent hadronic and 
semi-leptonic samples by first calculating the likelihood ratio estimator, Q = C{s + b)/C{b) using 
the likelihood functions from each method. We then create a combined estimator from the product 
of the semi-leptonic (Qsi) and hadronic (Qhad) likelihood ratio estimators, Q = Qsi x Qhad- The 
combined upper limit on the branching fraction is B{B^ — >■ t^i>t) < 3.3 x 10""^ at the 90% C.L. 

Using equation the combined branching fraction upper limit, and the measured value of \Vub\ 
m we set a limit on fs- We find fs < 0.480 GeV at 90% C.L. 
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Figure 6: The distribution of likelihood ratio as a function of B{B~^ — > T~^Ur). The dashed (solid) 
curves corresponds to the case in which the uncertainty on the expected background is included 
(not included). 

The BABAR Collaboration also searched for t^Vt decays using semi-exclusive semi- 

leptonic decay modes B~ — > D^i~v'iX where i = e, ^ and X can be a 7, vr*^, or nothing 9 . In 
the signal-side only the leptonic r decays, r"*" — > e~^Uel7r and r"*" — s- ^'^v^r, are identified. The 
B~ D*^i~V£ tags used in this current analysis overlap with the B~ — > D^i~i'^X tags and the 
analyses are therefore not combined. 

7 SUMMARY 

We have performed a search for the decay process B+ T^Vr- To accomplish this a sample 
of semi-leptonic B decays {D*^£^1ji) has been used to reconstruct one of the B mesons and the 
remaining information in the event is searched for evidence of B^ t^Ut- We find no evidence 
for this decay process and set a preliminary limit on its branching fraction of 

B{B+ T+Ur) < 4.3 X 10"^ at the 90% CL. 

By combining this analysis with a statistically independent B^ t^Vt search performed using a 
hadronic B reconstruction we find the preliminary combined limit: 

B{B+ T-^yr) combined < 3.3 X 10"^ at the 90% CL. 
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